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Abstract 

The work considers the damped Pmney equation, defined as the model 
arising when a hnear in velocity damping term is included in the Pin- 
ney equation. In the general case the resulting equation does not admit 
Lie point symmetries or is reducible to a simpler form by any obvious 
coordinate transformation. In this context the method of Kuzmak- 
Luke is applied to derive a perturbation solution, for weak damping 
and slow time-dependence of the frequency function. The perturba- 
tive and numerical solutions are shown to be in good agreement. The 
results are applied to examine the time-evolution of Gaussian shaped 
wave-functions in the Kostin formulation of dissipative quantum me- 
chanics. 

PACS numbers: 02.30.Hq, 02.30.Mv, 03.65.-w, 04.25.-g 

1 Introduction 

As is well-known, the Pinney equation [1] is ubiquitous in nonlinear dy- 
namics. A partial list of applications includes the exact solution for the 
classical and quantum harmonic oscillators [21 13], the search for invariants 
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(constants of motion) [U El E] , the stability analysis of charged particle mo- 
tion in accelerators [3 [S], the propagation of gravitational waves [H], the 
amplitude-phase representation of quantum mechanics [10], the derivation 
of the Feynman propagator for variable- mass problems [11], numerical solu- 
tions for non relativistic quantum problems [12], cosmological particle-crea- 
tion models [TB], cosmological models for the Friedmann- Robertson- Walker 
metric [2], isotropic, four-dimensional cosmological theories [13 [IS], rota- 
ting shallow water-wave systems [T7], curve flows in affine geometries [18], 
the stabilizer set of Virasoro orbits [19j, Bose-Einstein condensates with time- 
dependent traps and/or time- dependent scattering length [201 [21], discretized 
Pinney models [251 123] and nonlinear oscillations of transver sally isotropic 
hyperelastic tubes 

However, in spite of such a large number of applications, the Pinney 
equation in itself does not have any dissipation term. Hence it is natural 
to generalize the model by inclusion of the simplest damping mechanism, a 
term linear in the velocity. We call the resulting model the damped Pinney 
equation. As discussed in more detail in Section V, damped Pinney equa- 
tions arises in quantum mechanical models with dissipation. Specifically the 
Kostin formulation ^25j of the non conservative quantum time-dependent har- 
monic oscillator admits an exact Gaussian solution in terms of the solutions 
of a damped Pinney equation [26]. A damped Pinney equation also comes for 
Kostin's dissipative quantum mechanics under an arbitrary potential, after 
an expansion in the neighborhood of the classical path [27]. In addition, the 
Pauli equation for the Aharonov-Bohm effect with a time-dependent mass 
particle also reduces to a damped Pinney equation [28] . 

There is a close relation between the solutions of the time-dependent har- 
monic oscillator and of the Pinney equation, see the nonlinear superposition 
law at Section II. At first sight one can expect some simple connection 
between the solutions of the damped harmonic oscillator and damped Pinney 
equations. However, as is shown in the following, this relationship, if existing, 
is non trivial. Additional difficulties arises from the non existence of useful 
obvious point or nonlocal transformations removing the damping term and 
casting the damped Pinney equation into some well-known integrable class 
of equations. Reduction to Abel or Emden-Fowler equations can be shown to 
be trivial, but this does not help very much in solving the original problem. 
Indeed such classes of equations can be integrated only in some particular 
cases but not for the general damped Pinney equation. Finally the symmetry 
structure of the damped Pinney equation is poor in comparison to the three- 
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parameter group, SL{2,R), always admitted by the usual Pinney equation. 
These issues are discussed in more detail in Section II. 

In view of the above difficulties the present work is dedicated to the 
modest task of providing an approximate solution for the damped Pinney 
equation, assuming weak damping. In addition a slow variation of the fre- 
quency function is allowed. The perturbation technique chosen is the method 
of Kuzmak-Luke [221 EOl EI] , which is appropriate to generate bounded solu- 
tions for strongly nonlinear oscillator equations. Afterwards the perturbative 
solution is applied to dissipative quantum mechanics in Kostin's version. A 
damped Pinney equation appears naturally in this formulation of dissipative 
quantum systems. 

This work is organized as follows. In Section II the damped Pinney equa- 
tion is presented and some of its basic properties are discussed. In Section 
III the Kuzmak-Luke method is applied to derive a perturbative solution, 
free from secular divergences. Section IV shows good agreement of the ap- 
proximate solution with some numerical examples. In Section V the general 
connexion between the damped Pinney equation and dissipative quantum 
mechanics is explored. The final Section is dedicated to the conclusions. 



2 The damped Pinney model 

As is well known [L\, the general solution for the Pinney equation 

k 

X + uj'^(t)x = ^ , (1) 

X'^ 

where lo = Lo{t) is a time-dependent frequency function and k a numerical 
constant, can be written as 

X = {cial + C2(jI + 203^1^2)'/' , (2) 

where ci, C2 and C3 are constants such that C1C2 — c| = A; and cri and a2 are 
solutions for the time-dependent harmonic oscillator equation, 

(Ji + w^(t)(Ti = 0, i = l,2, (3) 

with unit Wronskian, cri6"2 — o"20"i = 1- We restrict k to be positive to prevent 
"collapse into the origin" issues. Notice the differences to the quantum case, 
where the wave-function for the time-dependent singular harmonic oscillator 
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is regular provided k > —fi^/{Am?)^ restoring dimensional quantities [52] . 
However, in the formal classical limit, setting ^ = we strictly need /c > to 
avoid X collapsing to the origin. This can be verified e.g. from the solution 
shown in Eq. ([2]) or from the potential function associated to Eq. ([1]), 
namely V = uj^x'^/2 + k/{2x'^) which is not bounded from below if /c < 0. As 
a particular example for = 1 and initial conditions x(0) = 1, x(0) = 0, Eq. 
([2]) gives = cos^ t + k sin^ t. Therefore one would have x = at some time 
t = > for any A; < 0. 

Equation ([1]) does not include any mechanism for dissipation. In this 
context it is natural to add a term linear in the velocity, yielding the damped 
Pinney equation 

k 

x + 2ex + uj'^{t)x = (4) 
x-^ 

where e > is a constant positive parameter. Notice that a time-dependence 
of the damping coefficient could be easily removed by an appropriate change 
of coordinates. Hence the basic properties of the system are already displayed 
for a constant e. 

Unlike the undamped (e = 0) case the integration of Eq. (jlj) is a chal- 
lenge in general. Indeed Eq. (jll) does not possess an universal Lie point 
symmetry group for arbitrary frequency function. This is in contrast to the 
richer sl{2,R) algebra always admitted by the usual Pinney equation [33] 
for any choice of uj{t). It can be shown that the damped Pinney equation 
has geometric symmetries for particular functional dependencies of uj{t), but 
these special cases are outside the scope of the present work. Also the ap- 
plication of symmetry generators which are quadratic in the velocity does 
not produce new results. Actually a coordinate transformation put Eq. (jlj) 
into the form of a generalized Emden- Fowler equation of index —3. There 
is an extensive literature on the integrability of generalized Emden-Fowler 
equations of arbitrary index |M1 135| , but we believe that it is still interesting 
to investigate in more detail the index n = —3 case. 

There is no easily identifiable point transformation, so that the damped 
Pinney equation could be put into a simpler, always integrable form. As a 
tentative example, consider the quasi-invariance transform 

x = p{t)Q{T), T = T{t), (5) 

where Q is the new dependent variable and p, T are functions of time to be 
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determined according to convenience. The resulting equation is given by 



pf' ^ + {pf + 2pf + 2epf) ^ + (P + ^V) Q = ^s- (6) 

From Eq. ([2]), it is easy to verify that it is not possible simultaneously to 
eliminate the damping term and set the coefficient of the inverse cubic term 
to a constant. 

Indeed the damping term in Eq. ([6]) can be eliminated provided 

-2et 

T=-^, (7) 
P 

ignoring an irrelevant multiplicative constant which could be included. By 
use of Eqs. (jUHZl) we derive 

where W = W{t) is a function of time, defined in terms of the auxiliary 
equation 

p + 2ep + u^p = — . (9) 

One can choose W = constant or even W = 0, but in all cases Eq. (jS]) 
would be still explicitly time-dependent through the inverse cubic term, since 
e 7^ and k is assumed constant. The non autonomous character of Eq. ([8]) 
prevents solvability in general. The time-dependence should be written using 
t = t{T), obtained via Eq. ([7]). 

At least, setting W = 0, one can cast the problem in the form of a 
generalized Emden- Fowler equation of index —3, 

where p(T) = A;exp(4et). This is equivalent to a choice of gauge [31]- How- 
ever, from the practical point of view of computing the solution, this property 
is not very helpful. Indeed consider the simplest case of a constant frequency 
function u = uq. Since ci; = 0, it is trivial to solve Eq. ([9]) assuming W = 0. 
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Using Eq. ([71) and the particular solution p = exp(— et) cosfycuo ~ e^t], Eq. 
dHD becomes 

— exp ( tan-^ Ju^ - 7^ T 1 . (11) 




Equation (|TT1) does not fall into any of the known integrable Emden-Fowler 
equations [351 EB]. Evidently, trying with different particular solutions for 
Eq. dn]) does not improve the scenario. 

For time-dependent frequencies the same procedure applies, but the re- 
sulting Emden-Fowler equation could be even more complicated. In the 
undamped (e = 0) case, Eq. (fTTl) reduces to the Ermakov-Pinney [H H] equa- 
tion, the integrability of which is well-known. Also in general the non au- 
tonomous Emden-Fowler equation of index —3 does not possess the Painleve 
property. For instance Refs. [511 E5] consider only the situation in which p 
is a constant. In addition there is no quadratic constant of motion for Eq. 
(HI), except in the autonomous case [3^ . 

A detailed analysis shows that the damped Pinney equation does not 
matches the Tresse-Cartan conditions [23 IM] and hence is not linearizable 
through a general point transformation. Also generalized Sundman trans- 
formations Q = Q{x,t), dT = F{x,t) dt, FdQ/dx 7^ 0, can be shown to be 
useless. Indeed the damped Pinney equation does not matches the conditions 
estabhshed in [391 SO] in order to be reduced to the free particle or similar 
simple equations via generalized Sundman transformations. 

Another avenue could be the reduction to an Abel equation of the second 
kind through v = v{x) = x, so that 

V— = -2ev - to^x + — . (12) 
dx x-^ 

However, Eq. ([T^ does not correspond to any of the known integrable Abel 
equations [HI |121 113] . Also certain classes of nonlocal symmetries are not 
helpful, see Eq. (5.10) of Ref. [Hj. The difficulties are due to the singular 
~ x~^ term. 

Some Pinney equations with time-dependent damping and nonlinear terms 
have been considered in the literature, e.g. in the calculation of the geome- 
tric phases and angles of dynamical systems. Specifically [13] these modified 
Pinney equations are given by 

X X + CO X = — —, (13) 
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where m, which can be interpreted as a time-dependent mass, and uj are 
arbitrary functions of time and is a constant. Notice that in Eq. fll3p the 
form of the time-dependence allows the use of 

« = ^, r=/§, (14) 



m p 

where p is an arbitrary function of time, to convert Eq. (fT3ll to a Pinney 
equation in standard form. 



P 



2m Am? ' 



k 



P+i^' + — -irz:^)pQ = -s- (15) 



Therefore these systems are somewhat trivial. Notice the particular form of 
the nonlinear term, which needs to be proportional to so that the scaling 
f|T^ works. From a physical point of view, the most relevant case would 
be the one in which the nonlinear term is time- independent. In addition 
other, different classes of generalized Pinney models with nonlinear and time- 
dependent damping have been discussed elsewhere [IB]. The conclusion of 
the Section is that apparently damped Pinney equations in the form (jl]) does 
not possess a closed form solution valid for general ivit), unlike the standard 
Pinney equation. The same applies to the alternative form given by Eq. (ITO!) 
for arbitrary p{t). 



3 Approximate solution for weak damping and 
slow time-dependence 

The difficulties enumerated in the last Section suggest the use of perturba- 
tion theory, in the case of weak damping. In other words we restrict to a 
small damping coefficient, e. An approximate solution would be welcome for 
physical applications as shown in Section IV. In this manuscript this modest 
goal of deriving an approximate solution is achieved through the method of 
Kuzmak-Luke [291 EOl EI] , which is good enough for equations which remain 
nonlinear even when the perturbation parameter goes to zero. Indeed this is 
the case of Eq. (jl]) when e 0. In addition the Kuzmak-Luke technique ap- 
plies only for slowly varying non-autonomous systems. Therefore we consider 
general frequency functions of the form 

u = n{et) , (16) 
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so that the time-dependence is slow. In Eq. fll6p . Q is an arbitrary ana- 
lytic function of its argument. To avoid some singularities it is assumed that 
n > 0. Fortunately Eq. (jl]) in the non perturbed (e = 0) case has a peri- 
odic solution so that all conditions to apply the Kuzmak-Luke method are 
satisfied. Naive perturbation theory would lead to a divergent solution. 

Other methods of removal of the divergences in the perturbation series, 
like the Lindstedt-Poincare approach, could also be chosen. However, we 
verified that the insistence on the use of such variational methods just adds 
unnecessary technical problems, even observing that Eq. (jl]) admits the 
Lagrange function 



L = L(x, X, t) 



^ I - 2 2 2 

[X — UJ X - 

2 ^ x2 



(17) 



In practice it can be checked that the Kuzmak-Luke method is among the 
simplest reliable approaches for the problem. 

The Kuzmak-Luke procedure seeks a series solution 



X 



Xo(r, t) + exi(r, t) + . . . 



(18) 



where t = et and 



dr 
'di 



fit) 



(19) 



with f = f{t) determined by the requirement that Xq be periodic in r with a 
constant period which can be taken as 27r without loss of generality. In the 
calculations r and t are regarded as independent variables. 
To zeroth-order in e, 



k 

Xq 



(20) 



with no damping term. Equation (120|) can be integrated once so that we 
obtain the slowly varying energy function 




H H 

^ 2 ^2x1 



The energy integral can be used to derive the quadrature 



1/2 



Xo = ± 



A\t) cos 



(21) 



1/2 



(22) 
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where the slowly varying amplitude A = A{t) and phase = are chosen 
according to convenience. For simplicity we ignore the time-dependence of 
the phase, setting (j) = constant In terms of the amplitude the energy is 
£^0 = f2vn^4M-fc, reproducing the usual law Eq ~ A'^ in the linear (/c = 0) 
case. From now we adopt the positive sign in Eq. fl22|) . 

By inspection, the trajectories defined by Eq. (l22l) become periodic with 
period 27r, independent on i, if we define 

f{i) = 2n{i), (23) 

We should determine A'^{t) requiring the first-order correction xi to be 
a periodic function of r, free from mixed secular terms. In the present case 
this can be shown to be equivalent to 

exp(2t) /-s- , / n\i)xl(T,i) k \ 

' ' drlEoit) ^ ^ °^ ' ^ — ^ = constant , (24) 



n{t) Jo V 2 2xl{T,t)^ 

where the integration is performed for fixed t. For more details see for in- 
stance Eq. (3.6.35a) of Ref. [29]. Eq. (IMD gives 

^2 _ Vn~oAoeM-i) ^ rV4|exp(-^ j ^^5) 

where Qq = ^(0) and Aq is a reference value. 

Finally the combination of Eqs. ( l22i) . ( l23i) and ( 125|) gives the zeroth-order 
solution 



1 

Xo = 



/r](et) 

/ Q /l2 -2et\l/2 . ( 

+ V^Aoe-^* j cos (^2 n{et')dt' 



1/2 



containing the two constants of integration Aq > and to. The perturbation 
procedure could be carried out to higher orders, but Eq. ( l26l) is sufficient for 
our purposes. 

The relevance of equation fl26l) cannot be overestimated. It provides a sort 
of JWKB solution for the damped Pinney equation, reproducing the quali- 
tative properties which are expected, namely, it shows a periodic motion 
spiraling toward the slowly varying "fixed" point k^/^/ n'^/'^{et). In addition 
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it can be checked that Xq > for all times provided k > 0. When there is no 
damping (e 0) nor nonlinearity {k 0), we derive Xq — > Aq cos[Qo{t — to)] 
as it should be. Of course one can adopt the simplistic viewpoint that a nu- 
merical solution could be sufficient so that the expression fl26|) is irrelevant. 
This naive criticism does not appreciate the usefulness of exact or approxi- 
mate analytical solutions in general. However, it is clear that expressions like 
d^n]) provide information not so easily found from numerics, as, for instance, 
the functional dependence of the solutions on the amplitude Aq. 

If there were no damping, but the frequency could still depends on i, Eq. 
( |26|) would be replaced by 



1 

Xo = 



n{et) 



^ + ^ + y^Ao (Vk + ^) 'cos (2^* n{et') dt'^ 



1/2 



(27) 

This is the form arising from the well-known exact solution ([2]) of the standard 
Pinney equation, considering the JWKB solutions ai = fi"^/^ cos[/ dt'] and 

= ^^^^'^ sin[/ ^Idt'] for the time-dependent harmonic oscillator along with 
appropriate parameters Cj. The powerfulness of both Eqs. fl26l) and fl271) is in 
their generality: they provide approximate solutions for the damped or stan- 
dard Pinney equations, irrespective of the functional form of the frequency, 
as long as it is slowly varying and positive. Near singular times, at which 
f2 = 0, the solutions would clearly be not appropriate. 

It is interesting to observe that 

i?o = l^(^^^^^ + v^) , (28) 

as follows from Eqs. fl?I]) and fl2B]) . generalizing the usual adiabatic theorem 
to the damped and nonlinear case. It is also consistent with the asymptotic 
approaching toward the "fixed" point /c^/"^/ where Eq = Vt\/k. 



4 Simple examples 

We have compared the approximate solution (126!) with numerical simulations 
for specific frequency functions as follows. 
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Figure 1: Approximate solution ( l26l) for i7 = fio and parameters e = 0.1, fio = 
k = 1,Aq = 2 and to = 0. Accordingly, the initial conditions are xo(0) = 
2.41,i;o(0) = —0.17. Direct numerical simulation reproduces the same curve. 



4.1 n = no 

For a constant frequency, Q = Qq, it is expected that the trajectory ap- 
proaches the value k^^"^/ fll^"^. Figure 1 shows the perturbative solution ( l26i) 
for the parameters e = 0.1,Qo = k = 1, Aq = 2 and to = 0. For the parame- 
ters, one found xo(0) = 2.41,i;o(0) = —0.17. Simulation of the corresponding 
damped Pinney equation with these initial conditions fully confirms the ac- 
curacy of the approximate solution with no noticeable difference to Fig. 1 . 

4.2 n = no {I + Py^/^ 

For a decaying frequency, = Qo{l + P)^^/^, the trajectory shows the 
asymptotic behavior xq ^^^A;^/^(l H-P)^/"^ because the confining po- 

tential becomes weaker. Figure 2 shows the perturbative solution ( !26l) for 
the same parameters as in Fig. 1. There is no need to show the numerical 
simulation, because it gives the same result as Fig. 2. 



5 10 15 20 25 30 

t 



Figure 2: Approximate solution fl26|) for 1] = fio (l + e^t^)~^/^ and parameters 
e = 0.1, f2o = k = 1,Aq = 2 and to = 0. Accordingly the initial conditions 
are Xo(0) = 2.41, xo(0) = —0.17. Direct numerical simulation reproduces 
the same curve. The dashed curve shows the asymptotic behavior x ~ (1 + 
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Figure 3: Approximate solution ( l26i) for f2 = f^o (1 + e^t^)^/^ and the same 
parameters of Figs. 1 and 2. The graphic is in accordance with the numerical 
results. The dashed curve shows the asymptotic behavior a; ~ (1 + e^t^)"^/'^. 

4.3 n = no{i + py/^ 

When the frequency grows as i7 = fig (1 + Y^'^i the trajectory shows the 
asymptotic behavior xq Q,Q^^'^k'^^'^{l+i'^ because the confining poten- 
tial becomes stronger. This is shown in Figure 3 with the same parameters 
as before. It reproduces very well the numerical simulation. 

To conclude the Section we have verified that the perturbation and nu- 
merical solutions disagree in the case of strong damping as can be expected. 

5 Application to dissipative quantum mechan- 
ics 

A popular approach for the time- dependent dissipative quantum harmonic 
oscillator is given by Kostin's model [25], which is expressed in terms of the 
Kostin equation 

2mdq^ \ 2 i \ip* ) ) dt ^ ' 
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n 



In Eq. fl29l) e is the damping coefficient and the remaining symbols have 
their usual meaning. In particular the wave- function is ip = %l){q,t). The 
Kostin nonlinear modification of the Schrodinger equation is among the only 
reliable alternatives to include damping in quantum mechanics, due to its 
smooth classical limit |17]. In particular it is better than considering non- 
Hermitian terms in the Hamiltonian since a continuity equation follows from 
Eq. m- 

Consider the de Broglie-Bohm decomposition ip = \Jn{q, t) exi){iS{q, t)/h), 
assuming the Gaussian Ansatz 

(..^W)-"%xp(-(i^)% (30) 

where x = x{t) and q = qci{t) are functions to be determined and the 
quantum-mechanical fiuid velocity is 

uiq,t) = = -{q- q^i) + q^i. (31) 
m oq X 

From the Kostin equation it follows [26] that 

x + 2ex + uj\t)x = (32) 
(id + 2eqci + uj^(t)qci = . (33) 

Clearly x satisfy a damped Pinney equation and qd solves the classical New- 
ton equation. The same equations can be derived for arbitrary nonlinear 
potentials, expanding around the classical trajectory [27]. Assuming weak 
damping and slowly varying frequencies, we can apply the results of Section 
III. 

As an example consider the oscillating frequency 

= fio(l +7sin(2et)) , (34) 

where f2o > and < 7 < 1 are fixed parameters. Using the approximate 
solution (12^ and numerically solving Eq. (155]) with the initial condition 
'?d(0) = l,gc/(0) = 0, we derive both the quantum fiuid density and velocity 
fields from Eqs. fl30]) and fl3Tl) . The results are shown in Figs. 4 (the 
standard deviation of the wave-packet), 5 (the particle density) and 6 (the 
velocity field at fixed position q = 0). In these graphs the parameters are 
7 = 0.7, h = m = Qq = l,e = 0.1,^0 = 4,to = 0. Other examples can be 
easily constructed as well. 
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Figure 4: Approximate solution for ( l32l) with uj = fio(l +7sin(2et)) and the 
parameters 7 = 0.7, h = m = Qq = l,e = 0.1, Aq = 4:,to = 0. 

6 Conclusion 

This manuscript derives an approximate solution for the autonomous damped 
Pinney equation for weak damping and slowly varying frequency using the 
Kuzmak-Luke method. The JWKB-like solution (!26l) reproduces with ac- 
curacy the numerical solutions of the model. Due to the usefulness of the 
Pinney equation in many areas of physics, mathematics and engineering, 
it would be relevant to derive general statements about the corresponding 
model when dissipation is present. In this context the perturbation approach 
here is just a first essay. However, probably there is no closed-form solution 
for the damped Pinney equation, valid for arbitrary frequency functions, in 
contrast to the undamped case. Adopting the optimistic view that such an 
universal solution could be constructed, one can expect the need for more 
complex nonlocal mappings than the generalized Sundman transformations. 
As a by-product there would be further insight on fundamental questions 
about dissipative non autonomous quantum mechanics. Up to now only 
special classes of damped Pinney equations or the equivalent Emden-Fowler 
equation ( JTOl) are known to be amenable to reduction of order. The interest- 
ing point about the expression (l26ll is that the frequency function is not of 
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Figure 5: The particle density as follows from Eq. (!30l) and the use of the 
approximate solution shown in Figure 4. Initial conditions for the auxiliary 
field: qciiO) = l,gd(0) = 0. 
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Figure 6: The velocity field associated to the same parameters of Figs. 4 and 
5, evaluated at q = 0. 

a particular form, except for having a slow time-dependence and not becom- 
ing zero. More accurate results could be found extending the perturbation 
theory to higher orders. 
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